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We consider how well the dark energy equation of state w as a function of red shift z will be 
measured using current and anticipated experiments. We use a procedure which takes fair account 
of the uncertainties in the functional dependence of w on z, as well as the parameter degeneracies, 
and avoids the use of strong prior constraints. We apply the procedure to current data from WMAP, 
SDSS, and the supernova searches, and obtain results that are consistent with other analyses using 
different combinations of data sets. The effects of systematic experimental errors and variations 
in the analysis technique are discussed. Next, we use the same procedure to forecast the dark 
energy constraints achieveable by the end of the decade, assuming 8 years of WMAP data and 
realistic projections for ground-based measurements of supernovae and weak lensing. We find the 
2a constraints on the current value of w to be Awo(2<r) = 0.20, and on dw/dz (between z = 
and z = 1) to be Awi(2<r) = 0.37. Finally, we compare these limits to other projections in the 
literature. Most show only a modest improvement; others show a more substantial improvement, 
but there are serious concerns about systematics. The remaining uncertainty still allows a significant 
span of competing dark energy models. Most likely, new kinds of measurements, or experiments 
more sophisticated than those currently planned, are needed to reveal the true nature of dark energy. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



I. INTRODUCTION 

Several years ago, observations of supernovae of type la 
(SNe la) demonstrated that the expansion of the universe 
is accelerating Associated with this acceleration is 

dark energy, a component of the universe with negative 
pressure, that makes up a significant fraction of the total 
energy density. Recent years have seen the supernova ev- 
idence for dark energy continue to mount 0, 0> IE IE 0> HI • 
Meanwhile, cosmic microwave background (CMB) data 
from the WMAP project [EDS, 

in combination with in- 
formation about either the Hubble constant ^lj or the 
galaxy power spectrum |l2t IT^ . also requires the exis- 
tence of dark energy. Cross-correlations between CMB 
anisotropics and matter power spectrum inhomogeneitics 
provide evidence for a late-time integrated Sachs- Wolfe 

(ISW) effect 0E EE1 BEI. This ISW effect id- 
eates a recent change in the inhomogeneous gravitational 
potential, providing further evidence for dark energy. 

One way to characterize the dark energy is to measure 
its equation of state, the ratio w = P/p of its pressure 
P and energy density p. The simplest model of the dark 
energy is a cosmological constant A, which has a constant 
w = — 1. Quintessence models describe the dark energy 
as a dynamical scalar field with an equation of state w > 
— 1, and in general, w will not be constant in time poll2ll 
I2E |23| . Furthermore, models of extended quintessence 
or "phantom energy" even allow w < — 1 |2J]. Models 
of the universe incorporating a cosmological constant or 
quintessence, in addition to cold dark matter, are known 
as ACDM or QCDM models, respectively. 

Discovering whether or not the dark energy is a cos- 
mological constant is the primary goal of the study of 
dark energy. If dark energy is shown not to be a cos- 



mological constant, then the next important issues are 
whether or not w < — 1, which can be theoretically prob- 
lematic HU (however, see [HH EHUHl for an alter- 
native viewpoint), and whether or not w is changing with 
time. While a quintessence with constant w ~ — 1 is very 
difficult to distinguish quantitatively from a cosmologi- 
cal constant, the qualitative difference for fundamental 
physics is enormous, so it is critical that the maximum 
effort be made to reduce the uncertainty in | w + 1 1 and its 
time derivative. Ultimately, the precise quantitative val- 
ues of w and its time derivative are important for model 
building, but this is less crucial, at present, than the qual- 
itative issue of whether the dark energy is dynamical or 
not. In order to address such questions about the nature 
of the dark energy, it is crucial that the equation of state 
and its time variation be determined observationally. 

It has been known for some time that the cosmological 
probes used for stuyding the dark energy are plagued by 
numerous parameter degeneracies [3l|, |32|, 13a . 134 . 1351 l3a . 
l37| . Because of these degeneracies, a large uncertainty 
in a cosmological parameter not directly related to the 
dark energy sector (e.g., the matter density) can lead to 
a large uncertainty in the dark energy equation of state. 
This sensitivity to other parameters means that one must 
incorporate the uncertainties in all parameters to obtain 
a realistic estimate of the uncertainties in the dark energy 
equation of state. In particular, one should avoid the use 
of strong priors on the form of w(z), the values of other 
parameters, and independent information coming from 
other experiments. These can lead to underestimates of 
the uncertainty in w by a very large factor. We will also 
discuss how standard likelihood marginalization can give 
a misleading impression of the uncertainty in w. 

The analysis presented here employs a x 2 minimiza- 
tion procedure in order to avoid possible problems 
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with marginalization, and assumes only weak prior con- 
straints. We compare our x 2 minimization with the stan- 
dard approach, the marginalization of a probability func- 
tion computed using Monte Carlo Markov Chains. It is 
argued that our procedure gives a more conservative as- 
sessment of constraints in certain degenerate parameter 
spaces. In this sense, our procedure and the Markov 
Chain are complementary, as we will discuss. 

We apply this \ 2 minimization analysis to the cur- 
rently available data. The simultaneous determination of 
the dark energy equation of state and its red shift deriva- 
tive is difficult without the combination of data from sev- 
eral cosmological probes. Thus we analyze the latest data 
from WMAP, SN la searches, and the Sloan Digital Sky 
Survey. Strong priors are avoided to prevent the under- 
estimation of dark energy constraints. Although the data 
favor w(0) < —1 and ti/(0) > 0, we are unable to rule out 
the cosmological constant conclusively. Furthermore, we 
find the constraints on w(z) to depend strongly on the 
parameterization chosen for w(z), as well as on the loca- 
tion of the best-fit model in parameter space, despite the 
fact that we have restricted our study to two-parameter 
equations of state. We also list possible biases in the 
cosmological probes, and discuss their effects on our re- 
sults. Our conclusion that the cosmological constant is 
consistent with current data agrees with previous analy- 
ses, as we discuss in Sec. II VI This establishes the validity 
of our analytic methods, which we can then apply with 
confidence to projections of future measurements. 

Since the cosmological constant is consistent with cur- 
rent data, it is useful to ask precisely how well dark en- 
ergy may be constrained in the future using known cos- 
mological probes. We address this question by applying 
our x 2 minimization procedure to simulations of the data 
from these probes. Besides the CMB and SNe la, we 
simulate data from a probe of weak gravitational lensing 
(WL). In the future, WL is expected to be useful for con- 
strainingdark energy ©Ellllllllllllllllllllli 
l47l Esl US IHoL l5l| . Most of these studies have demon- 
strated that WL breaks parameter degeneracies in the 
CMB analysis, and have shown the effectiveness of the 
combination of CMB and WL data in constraining dark 
energy. 

Our study is the first to use a joint analysis of simu- 
lated CMB, SN la, and WL data, rather than imposing 
prior constraints as a "replacement" for one of these data 
sets. This is especially important for parameters such 
as w;(0) and w'(0), which have strong degeneracies with 
other parameters. Since our goal is the study of w{z) 
constraints, rather than general parameter constraints, 
we assume the cosmological model outside of the dark 
energy sector to be as simple as possible. Within this 
model, we are careful to avoid strong priors on the cos- 
mological parameters. It is well understood that each of 
these cosmological probes has systematic uncertainties 
associated with it, and that significant progress needs to 
be made in order for each one to reach the level of pre- 
cision necessary for constraining dark energy. Therefore, 



we attempt at every opportunity to strike a reasonable 
balance between optimism and realism. 

The paper is organized as follows. Sec. Ill Al discusses 
several parameterizations of the dark energy equation of 
state. Our choice of cosmological parameters, and the 
prior constraints imposed on them, are listed in Sec. Ill Bl 
X 2 minimization is compared with marginalization in Sec. 
Ill CI Sec. IIII Al summarizes our analysis of current data, 
and Sees. IIII BlIIlTjl describe the simulation and analysis 
of future data. Finally, Sec. IIVI lists and interprets our 
findings, for our analysis of current data as well as our 
forecasts of future constraints, and Sec. discusses our 
conclusions. 



II. ANALYSIS METHODS 
A. Dark energy parameterization 

In the spirit of reasonable optimism, we have as- 
sumed that the dark energy may be parameterized us- 
ing an equation of state w(z) = P/p, where P and p 
are the pressure and dark energy density, respectively, 
and w{z) is an unknown function of red shift. We 
have also assumed that the dark energy sound speed is 
c s 2 = 1. Without any theoretical guidance about the 
form of w(z), the space of all possibilities is equivalent 
to an uncountably infinite set of cosmological parame- 
ters. The data depend on w(z) only through a multi- 
ple integral relation |3?1 |. so we lack the large number 
of measurements of w(z) that would be necessary for 
non-parametric inference |H^ . Thus we must describe 
w(z) using a small number of parameters in order to 
keep the analysis tractable. The danger is that there 
is no "natural" way to parameterize w(z), so the choice 
of parameterizations is essentially arbitrary. Ideally, one 
would like to verify that any constraints obtained on w(z) 
are parameterization- independent. However, the space of 
all possible parameterizations is infinite-dimensional, so 
we're back to square one. The only feasable option is 
to test for parameterization-independence by comparing 
constraints on a small number of "reasonable-looking" 
parameterizations. 

To begin with, we assume that w{z) is an analytic func- 
tion, and that it can be approximated at z <C 1 by the 
first few terms in its Taylor series about z = 0. All 
known probes are sensitive to an integral of some func- 
tion of w(z), rather than to w(z) itself, so any "bumps 
and wiggles" in w(z) are effectively smoothed out. It 
is reasonable to assume that this smoothed w(z) can be 
approximated by a low-degree polynomial in z out to 
2 <i 1. From now on we proceed based on this assump- 
tion, with the caveat that the "true" equation of state can 
have added "bumps and wiggles" that disappear when 
smoothed. (Note, however, that any large-scale kinks 
or sharp changes in the actual equation of state can be 
missed by such simple parameterizations. For an alter- 
native approach, see, e.g., H3-) 
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TABLE I: Summary of the dark energy equation of state pa- 
rameterizations used in the literature. 



The simplest parameterizations useful for describing 
dynamical dark energy have at least two parameters, 
one describing the equation of state w(Q) at z = and 
the other parameterizing its red shift derivative w'(0) = 
Q . A few parameterizations in the literature are 



d 

dz 



w(z) — wq + wz (Simple parameterization) (1) 
w(z) = w + (1 - a(z)) w a 

= w ° + TT~ ( Refs - SHI) (2) 
w(z) = w + (1 - a(z)) w a + (1 — a(z)) 2 w b 

WnZ WhZ 2 ,„ I 1 % , . 

wo + „ ° ; + ^ b „ (Ref. 56]). (3) 



(1 + z) (1 + z) 2 

(Note that w a in © differs by a factor of —1 from the 
definition introduced in j56j. We chose 10 in order to 
facilitate comparison with © and (0J.) We are interested 
in the simplest two-parameter equations of state, and as 
we will see, Q is of limited utility for studying high 
red shift data. This leaves only Eq. @ out of the three 
parameterizations above. Since we would like to compare 
two "reasonable" two-parameter equations of state, we 
introduce a fourth parameterization, 



w(z) 



wo + w\z if z < 1 
Wo + wi if z > 1. 



(4) 



Dark energy constraints based on (@J and (J2J) may be 
compared in order to estimate the parameterization- 
dependence of our results. Table [I] compares all four pa- 
rameterizations using their present-day derivatives u>'(0), 
their mean low red shift (0 < z < 1) derivatives w 1 = 
w(l) — w(0), and their high red shift limits lim 2 _ (00 w(z). 

The high red shift limit of each parameterization be- 
comes important when considering data from the CMB. 
Rather than providing information about w(z) directly, 
the CMB data require only that the dark energy den- 
sity be less than ~ 10% of the critical density at the 
red shift Zdec ~ 1100 of photon decoupling |53,|53|. If 
w(z) > 0, the dark energy density will increase faster 
with red shift than the matter density, making the CMB 
constraint extremely difficult to satisfy. For the four dark 
energy parameterizations discussed above, the CMB ef- 
fectively requires that w(z) < at red shifts z> 1, that 
is, that the quantities listed in the fourth column of Table 
[I] are no greater than zero. 
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FIG. 1: Upper bounds on w(z) for the three equations of state 
and (J2J. In all three cases, too has been fixed at —1. 



In each of the the two- and three-parameter equations 
of state QEIOIQJ, this high red shift constraint prevents 
w{z) from changing too rapidly at low red shifts. For 
fixed wo, we see that: w' < in Q; w\ < —wo in 
and w a < —wo in (J2J. In particular, w\ and w a have the 
same upper bound. Note that these bounds arise purely 
from the choice of parameterizations for w(z). There 
is no fundamental physical reason for assuming that an 
equation of state which becomes positive at red shift z = 
1 will stay that way long enough to interfere with the 
physics at decoupling. 

Fig. n snows the upper bounds on the function w(z) 
for each of the three parameterizations , J5J , and Q , 
assuming wo = —1 in each case. Parametrization 
disallows essentially all positive w'; we will not study 
it further. Meanwhile, J2J allows w(z) to increase by 
w' = w(l) — w(0) = 0.5 at low z, while allows w(z) to 
increase by w' = 1. Some parameter combinations such 
as (wo = —1, w' — 0.6) are acceptable in (@J, but are 
simply not allowed in @. That is, the two parameter- 
izations cover different regions in the (w, z) plane. Of 
course, these parameterization-dependent effects of the 
high red shift w(z) constraint will be important only if 
the dark energy is found to have w(z) ~ at z > 1. 

We choose to use parameterization Q for most of the 
subsequent work, since it allows the widest range of w' 
values. Equation J5J is also studied, and the results based 
on the two parameterizations are compared in order to 
estimate the parameterization dependence of our dark 
energy constraints. We describe the dynamics of the dark 
energy in terms of the ratio of the dark energy density to 
its value today, 



Q(z) = Pde {z)/p de (0). 



(5) 



For parameterization |0J, the fluid continuity equation 
can be used to show that 



(1 + 2)3(1+^0-^1) e 3w lZ if z <l 5 
(1 + z )3(l+«>o+t«l) e 3u!i(l-21n2) Jf ^ > J. 



(6) 
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(A corresponding: exp ression for equation of state 10) can 
be found in j5J, log.') Given our assumption of a flat 
universe, the Hubble parameter H{z) = a/a evolves with 
red shift as follows. 

H(z) = H oy /Sl m (l + z)3 + (l-n m )Q(z) (7) 

These results will be used in our discussions of the cos- 
mological probes. 

B. Cosmological parameters, degeneracies, and 
priors 

Cosmology may be described using a large number of 
parameters corresponding to a wide range of possible ef- 
fects. The approach taken by, e.g., 0] is to test for as 
many of these effects as possible; they describe cosmology 
in terms of 13 parameters. On the other hand, our aim is 
more specific. We wish to determine whether dynamical 
dark energy will be distinguishable from a cosmological 
constant by the end of the decade. If we can answer this 
question in the negative after considering only a subset 
of these 13 parameters, then the addition of more pa- 
rameters will not change our conclusions. We find that 
even with a relatively limited set of parameters we will 
not be able to rule out dynamical dark energy. Thus, in 
the interests of simplicity and computational efficiency, 
we restrict ourselves to nine parameters. 

Outside of the dark energy sector, we choose the sim- 
plest possible description of the universe that is con- 
sistent with observations. The universe is assumed to 
be flat, and we do not include tensor modes, massive 
neutrinos, or primordial isocurvature perturbations. Be- 
sides the dark energy parameters Wq and W\ , our cosmo- 
logical models are parameterized using: h = Hq/(100 
km sec -1 Mpc" 1 ), where Hq is the Hubble constant; 

LO m = Q m h = p m h I p ar it\ U>b = &bh — Pbh I Peril, T, 

the optical depth to reionization; A, the normalization of 
the CMB power spectrum; n S} the spectral index of the 
primordial power spectrum; z s , the characteristic weak 
lensing source red shift. (The WL source red shift was 
shown to be important in, e.g., |59l l6Ct |6l|.l Actually, 
[l3| shows that n s is not absolutely necessary, since the 
simple Harrison-Peebles- Zeldovich spectrum (n s = 1) fits 
the data well. However, n s is a well-motivated parame- 
ter whose degeneracies with the dark energy parameters 
could be important. Borrowing the terminology of |13| . 
we call (h, uj m , ojb, t, n s , A) the "vanilla" parameters. 
Thus, our analysis adds to the simple "vanilla" model 
one parameter describing weak lensing sources, and two 
parameters describing the dark energy. 

Our constraints based on this limited set of parameters 
will be optimistic, since including more parameters will 
tend to increase the uncertainties in wq and W\. This 
is especially true when the parameters that we neglect 
are highly correlated with some of the parameters that 
we do consider. Ref. p| shows a strong degeneracy 
between the Hubble parameter h and the curvature 



that could potentially affect our dark energy constraints. 
Also, 39] points out that neutrinos become important 
when using weak lensing in the nonlinear regime to study 
dark energy. 

h(w ,w 1 ) 



1 r 




FIG. 2: h(u)Q,W\) satisfying the CMB angular diameter dis- 
tance degeneracy, with parameters u m = 0.14, Zdec = 1089, 
and d% c) = 14.0 Gpc fixed to their WMAP values. 



Even in our relatively simple parameter space, there 
exist several parameter degeneracies. That is, to each 
point in parameter space there corresponds a degenerate 
region of observationally indistinguishable points. A de- 
generacy that is particularly important to the study of 
dark energy is the angular diameter distance degeneracy 
H3- E3 US EH • For flat universes containing matter and 
a dark energy of the form this implies that models 
with the same values of ui m , ^b, and the comoving angu- 
lar diameter distance d^cf 1 to the decoupling surface, will 

produce indistinguishable CMB power spectra, as 
a function of h, wq, and w\ is given by 

(dec) r d " dz' 

d AC (h,w ,wi) = cj W^rj' ( 8 ) 

where H(z') is given by l(7jl. Even when u) m is fixed, each 
choice of wo and wi has a corresponding h for which this 
degeneracy is satisfied, as shown in Fig. [21 Other im- 
portant degeneracies include the large-scale CMB degen- 
eracy |3l[ among oj m , t, n s , and w(z); the small-scale 
CMB degeneracy [33| among t, A, andn,; and the su- 
pernova luminosity distance degeneracy [33 between O m 
and w(z). The presence of these parameter degenera- 
cies means that we must combine several cosmological 
probes in order to obtain reliable constraints. Further- 
more, combinations of known probes contain residual de- 
generacies, which we must deal with carefully in order to 
determine constraints on the dark energy. 

In order to make a fair assessment of parameter degen- 
eracies, it is necessary to avoid strong prior constraints 
on the cosmological parameters. We have been careful to 
use very weak priors. The dark energy equation of state is 
required to be negative or zero at all red shifts, w(z) < 0. 
For the model this implies the two prior constraints 
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WO + Wl 


N/A 






TABLE II: Summary of the prior constraints assumed for the 
X 2 minimization analysis. 



wo < and wq + wi < 0. For physical reasons, we must 
assume < Cl m < 1, u> b < u) m , r > 0, and z s > 0. Be- 
sides these, the constraints 0.4 < h < 1.1, 0.003 < uj b , 
uj m < 1, r < 1, 0.5 < A < 1.5, 0.5 < n s < 1.5, and 
z s < 1.5 are imposed for computational convenience. The 
above priors are summarized in Tabic ITT1 

The simulated data sets used here for constraint fore- 
casts assumed a fiducial model {wq = —1, Wi = 0, 
h = 0.7, oj m = 0.15, oj b = 0.023, r = 0.1, A = 0.8, 
n s = 1, z s = 1}. This roughly corresponds to the "vanilla 
lite" ACDM model of [ij, with lensing sources assumed 
to be at characteristic red shift 1. Our chosen fiducial 
model implies f2 m = 0.31, f2b = 0.047, and as = 0.92. 

C. x 2 minimization 

In order to provide an accurate picture of the dark 
energy constraints in a degenerate parameter space, a \ 2 
minimization procedure was used. Constraints on the 
dark energy parameters were determined using the la 
and 2a contours of the % 2 function in the (wq, wi) plane. 
At any given point (wq, wi), x 2 (wo,wi) was computed 
by minimizing over all other cosmological parameters. 

X 2 minimization handles certain types of degenera- 
cies more carefully than marginalization, the proce- 
dure most commonly used to obtain cosmological con- 
straints, meaning that the two methods are complemen- 
tary. Marginalization favors models that fit well over 
"large" regions of parameter space. Given two cosmolog- 
ical models with the same low x 2 , marginalization will 
give extra weight to the model that sits in a large region 
of low-x 2 models. Meanwhile, % 2 minimization will pick 
the best model while completely ignoring its surround- 
ings. These differences are illustrated in Fig. |3 based on 
a sample x 2 function of two arbitrary parameters x and 
y. The la contour of X 2 ( x tV) is shown in Fig. E^top). 
Note that the contour has a small slope at low x, but is 
nearly parallel to the y axis at high x. 

If we are interested only in constraints on x, then we 
can either minimize or marginalize over y. Fig. I^Jmid) 




Bounds on x: 

marginalized: 
minimized: 



FIG. 3: Comparison of minimization and marginalization in 
a (hypothetical) degenerate parameter space. Two models, 
A and B, are labelled, (top): Sample la contour of x 2 ( x yV) 
in a degenerate parameter space, (mid): Probability func- 
tions based on minimization and marginalization. The prob- 
ability based on minimization (dashed line) is V m in(x) oc 
exp (mhiy x 2 ( x i y))- The marginalized probability (solid line) 
is P marg (x) oc J exp ( y x 2 (x,y))V pr ior(x,y)dy. (hot): Con- 
straints on x based on minimized and marginalized probabil- 
ity functions. 

sketches the probability functions obtained using the two 
different methods. Consider the two models A and B, at 
xa and xb respectively, whose y values are chosen to 
minimize x 2 (x,y) at the corresponding x values. Mini- 
mization takes the viewpoint that, since \\ ~ Xb> the 
x values xa and xb are approximately equally proba- 
ble. On the other hand, marginalization assigns a much 
greater probability to xb, since there are many more low- 
X models at xb than at xa- Thus the x 2_mm imizcd 
probability distribution looks like a broad plateau, while 
the marginalized probability looks like a sharp peak with 
a "shoulder." The bounds on x obtained from these two 
methods, as shown in Fig. EHbot), are quite different. 
Note that this discussion is not just academic. The CMB 
95% probability contour in the (h, fltot) plane, shown in 
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TABLE III: Comparison between na \ 2 contours and proba- 
bility contours for a Gaussian likelihood function of two vari- 
ables. 

Fig. 7 of ^3 ' i s qualitatively very similar to the contour 
shown here in Fig. Bltop). 

The parameter bounds derived from marginalization 
have excluded models such as A, which are within the la 
X 2 contour, simply because the low-% 2 regions at their x 
values are "thin" in the y direction. A procedure that 
excludes such "thin-contour" models, in favor of "thick- 
contour" models such as B, can be problematic for two 
reasons. First of all, it is not clear that low-% 2 mod- 
els such as A really should be ignored; it is prudent, at 
least, to know of the existence of such a model. Points 
at the ends of long, narrow regions of degeneracy are 
"thin-contour" points, so marginalization can underesti- 
mate the sizes of these degenerate regions. Such points 
should be considered in a proper treatment of parameter 
degeneracies. 

Secondly, the notions of "thinness" , "thickness" , and 
"distance" in parameter space depend on the prior con- 
straints chosen. For a parameter vector p of N param- 
eters, the a posteriori probability marginalized over all 
but m of the parameters is given by 

p(p°, . . . ,p m -') = J d P m . . . d P N - 1 c(p)r p „ or (p). (9) 

It is customary to reduce this dependence on priors by 
choosing uniform priors, with V P rior{p) constant on a 
subset of parameter space and zero elsewhere. However, 
the concept of uniform priors is itself parameterization- 
dependent. Choosing priors that are uniform in a differ- 
ent set of parameters q is equivalent to letting V pr ior{p) 
equal the Jacobian determinant of the transformation 
from p to q. 

On the other hand, marginalization has several advan- 
tages of its own. First of all, one can always come up 
with a X 2 (. x >y) function which is handled correctly by 
marginalization, but for which minimization implies arti- 
ficially tight constraints. Also, marginalization naturally 
defines a probability in parameter space. An analysis 
based on marginalization is capable of providing proba- 
bility contours, which may be more useful than x 2 con- 
tours. If the likelihood function is Gaussian, these two 
types of contours are related in a simple way, as shown in 
Table ITTT1 However, for a non- Gaussian likelihood func- 
tion, there is no simple way to relate contours of X 2 (p) 
to a given probability contour without marginalization. 
Meanwhile, it may be the case that a "natural" set of 



parameters exists for describing a theory, such as h, f2 m , 
and Qb in general relativistic cosmology. The existence of 
"natural" parameters makes it easier to choose a "reason- 
able" prior probability distribution V pr ior, since one no 
longer needs to worry about reparameterization changing 
the form of Vprior- Constraints derived from a marginal- 
ization over such parameters, assuming weak and uniform 
priors, can be quite convincing. 

We take the point of view that any final constraints on 
wo and w\ should be independent of the analysis method 
used. If a claim made using one of the two methods does 
not hold up to scrutiny by the other method, then the 
issue is too close to call. For example, we do not be- 
lieve that model A in Fig. [3|is ruled out, even though 
it is excluded by marginalization over y. In this sense, 
X 2 minimization and marginalization are complementary 
analysis techniques that are useful for handling different 
types of parameter degeneracies in a non-Gaussian likeli- 
hood function. We have chosen to use x 2 minimization, 
partly because we are concerned about degenerate re- 
gions such as in|S{top), and partly because it is not clear 
which parameters should be used to describe dark en- 
ergy, reionization (r or z re ion), and the power spectrum 
amplitude {A or erg). 

For completeness, we discuss a few computational is- 
sues here. Minimizations were performed using the 
Amoeba routine of 62], although the Powell routine from 
that reference was found to give very similar results. 
The computation was sped up by generating CMB power 
spectra for 100 different n s values at once, and interpo- 
lating to find x 2 f° r intermediate n s values. This gave us 
information about x 2 over our entire range of n s values 
with just one call to CMBFAST. (Interpolation-related 
errors in x 2 were found to be negligible.) Also, we min- 
imized separately over the three parameters A, n s , and 
z s , since a variation in one of these parameters did not 
require the recomputation of CMB power spectra. This 
minimization over A, n s , and z s was nested within the 
minimization over h, uj m , <-^b, and r. Using this technique, 
Xcmb could be minimized at one point in the (wq, wi) 
plane in about 2 — 4 hours, on a single node of the Hydra 
computer cluster. 



III. COSMOLOGICAL PROBES 
A. Current data 

Since our analyses of current data follow standard pro- 
cedures, wc will discuss them only briefly. 

Supernovae of Type la: Supernovae of type la (SNe la) 
are standardizable candles jSJ l64l l65l |6£| , with magni- 
tude m(z) at red shift z given in, e.g., |67j, by 

m = 5 \og w {D L )+M. (10) 

The dimensionless luminosity distance Dl(z) = (1 + 
z)Hq Jq dz' I H{z') depends on f2 m and the dark energy 
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parameters through H (z) , given in (J7J) . Our analysis min- 
imizes x 2 with respect to the magnitude parameter Ai, 
which is dependent on Hq and the SN la absolute mag- 
nitude. 

Cosmic Microwave Background: For the dynamical 
models of dark energy considered here, power spectra 
were computed with a modified version of CMBFAST 4.1 
[68|. provided by R. Caldwell Our analysis of the 

data used our own implementation of the CMB \ 2 func- 
tion described in [jjj. Noise parameters and constants 
describing the WMAP parameterization of the Fisher 
matrices were taken from the data tables created by [tT| 
and and provided with the WMAP likelihood code. 

Galaxy Power Spectrum: Our analysis of Sloan Digital 
Sky Survey (SDSS) data used the galaxy power spectrum 
HH and the likelihood code of 0] , made publically avail- 
able at . Given a cosmological model, the correspond- 
ing matter power spectrum was computed using CMB- 
FAST. The normalization of the power spectrum was 
treated as a nuisance parameter; \ 2 was minimized with 
respect to it. An accurate formula for the matter power 
spectrum in the nonlinear regime, for dynamical dark 
energy cosmologies, is not currently available. Further- 
more, galaxy biasing in the nonlinear regime is even less 
well understood. Therefore, the SDSS data analysis pre- 
sented here only used measurements for k < 0.15/i/Mpc, 
as recommended in 131. 



B. SN la Simulation 

1. Supernova simulation strategies 

Monte Carlo techniques were used to simulate the mag- 
nitudes and red shifts of Type la Supernovae, from pre- 
vious, current, and future supernova surveys, that would 
be available for analysis by the end of the decade. Rather 
than mixing real and simulated data, it was decided to 
simulate all of the surveys from scratch. In order to con- 
duct a realistic simulation of future supernova surveys, 
we studied the available literature, including analyses of 
current data as well as plans for future surveys. For each 
SN la survey we estimated the expected number of su- 
pernovae to be observed and modeled a SN la red shift 
distribution. The final simulated SN la data set included 
2050 supernovae ranging in red shift from 0.01 to 2.0, 
as summarized in Table IIVI Our simulation was based 
upon the following assumptions, which we justified by 
reference to data whenever possible. 

• 25% of the SNe la found by each survey were as- 
sumed to be unusable for cosmological analysis. 
This is based on the fact that 23% of the SNe la 
in the Tonry/Barris data set 00 did not survive 
their cuts on very low red shifts (z < 0.01) and 
galactic host extinctions. 

• Two-thirds of all low red shift SNe (of all types) 



were assumed to be of type la, based on the super- 
novae reported in [t3 | . 

• Red shift distributions were assumed to be uniform 
for low red shifts. Since m{z) at low red shifts 
is independent of cosmological parameters other 
than M., the details of the red shift distribution 
for z < 0.1 should be unimportant for parameter 
constraints. 

• Red shift distributions were assumed to be Gaus- 
sian (Piz) oc exp ^— ) for higher red 

shift surveys, unless otherwise specified in the lit- 
erature. 

• Current data reported in IE were simulated 
by dividing them into three "surveys:" 

- low z, Tonry/Barris SNe la with z < 0.1; 

- mid z, SNe la from 00,0 with z > 0.1; 

- high z, HST GOODS SNe la 0. 

This was done because the current data come from 
a collection of previous surveys over a wide range of 
red shifts, and are poorly approximated by uniform 
or Gaussian distributions. 

• 40% of the expected Dark Energy Camera [75| data 
were assumed to be available, bringing the total 
number of supernovae to just over 2000. Even un- 
der optimistic assumptions about systematic uncer- 
tainties, we do not expect an increase in the number 
of SNe beyond this point to improve constraints on 
dark energy 

2. Magnitude uncertainties 

A well-studied SN la sample, with accurate spectro- 
scopic measurements of the supernova host galaxy red 
shifts, can have an intrinsic magnitude uncertainty for 

each supernova as low as dm ^ — o.i5 [zllz3- 

If less 

accurate host red shift information is available, the mag- 
nitude uncertainy suffers; |75| estimates an uncertainty 

(int) 

of 0Vn ~ 0.25 when only photometric red shift infor- 
mation is available for the host galaxies. Meanwhile, the 
average magnitude uncertainty in the Tonry/Barris sam- 
ple, for SNe la in the red shift range 0.1 < z < 0.8, is 
at™'' = 0.3. Since the final dataset will be a combi- 
nation of many supernovae with varying amounts of red 
shift information, the approach adopted here is to assume 
a magnitude uncertainty of (Tm"*' = 0.2 for supernovae 
in the intermediate red shift range 0.1 < z < 0.8. 

Judging from the supernovae tabulated in and 0, 
the low red shift supernovae tend to have smaller mag- 
nitude uncertainties than average. The mean magnitude 
uncertainty is 0.18 for the z < 0.1 SNe la, compared to 
0.25 for the full Tonry/Barris sample. Thus, for the pur- 
poses of this simulation, it is assumed that the magnitude 
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FIG. 4: Supernova apparent magnitude uncertainty used in 
the simulation 



uncertainty for each low red shift (z < 0.1) supernova is 

the minimum value, Um — 0.15. 

Meanwhile, supernovae at red shifts greater than about 
0.8 — 1.0 become increasingly difficult to observe from the 
ground, and the spectroscopic analyses of such SNe la be- 
come very time-consuming [Zllli In the Tonry/Barris 
sample, the magnitude uncertainty rises from a mean of 
0.3 in the red shift range 0.1 < z < 0.8 to 0.35 in the 
range z > 0.8. Based on this, it was decided to assume 



(int) 



an intrinsic magnitude uncertainty of a 
supernovae with z > 0.8. 

The correspondence (|10f) between apparent magnitude 
and red shift is exact only in the absence of a peculiar 
velocity. It is necessary to include an extra uncertainty 
to account for the nonzero peculiar velocities of actual 
supernovae. Following pj, it was assumed here that the 
peculiar velocity uncertainty of each supernova was a v = 
500 km/ sec. The fractional uncertainty in the luminosity 
distance due the peculiar velocity uncertainty was taken 
to equal the fractional uncertainty ^j^- in red shift. The 
resulting contribution to the magnitude uncertainty is 



(pec) _ 



ln(10) cV 



(11) 



where c is the speed of light. The resulting form of the 
magnitude uncertainty, after adding (Jm ec ^ in quadrature 
to the intrinsic uncertainty, is plotted in Fig. 0] 

Finally, following [7(| , we assumed a systematic uncer- 
tainty that prevents constraints on the apparent mag- 
nitude m(z) from becoming arbitrarily tight. In the 
absence of such a systematic effect, the magnitude un- 
certainty in a bin containing Nun supernovae will be 
a m/ \/~{Nbin), which approaches zero with increasing 
Nbi n . To this we add 5m = 0.04 in quadrature, which 
imposes a floor 5m on the uncertainty in each bin. We 
choose the bin size Az = 0.1; this acts as a "correlation 
length" for the systematic. For a single supernova in a 
bin containting Nbin supernovae, this systematic implies 



SN la 


number 


red shift 


data set 


or olNe la 


distribution 


Current (low z) 


83 


uniform, 






0.01 < z < 0.1 


Current (mid z) 


117 


Gaussian, 






Zmean = 0.55, O z = 0.25 


Current (high z) 
[8] 


16 


Gaussian, 

Zmean = 1.0, O z = 0.4 


Carnegie (low z) 


94 


uniform, 


[79] 




0.01 < z < 0.07 


Carnegie (high z) 


90 


Gaussian, 


[79] 




Zmean = 0.4, <7 Z = 0.2 


DE Camera 


570 


Gaussian, 


[75] 




Zmean = 0.55, O z = 0.25 


ESSENCE 


150 


uniform, 


[80, 81, 82] 




0.15 < z < 0.75 


PANS 


80 


Gaussian, 


[83, 84] 




Zmean = 1.0, CT Z = 0.4 


SDSS 


100 


uniform, 


f85] 




0.05 < z < 0.15 


SNfactory 


225 


uniform, 


[7J 




0.01 < z < 0.1 


SNLS 


525 


Gaussian, 


[77, 86, 87, 88] 




Zmean — 0.6, tj z — 0.3 


Total 


2050 





TABLE IV: Numbers and red shift distributions of simulated 
type la supernovae. 



25 for an effective magnitude uncertainty 



^ //) = y(^ nt) ) +(^ pec) ) +N Un 5m?, (12) 
which increases with increasing Nbin- 

3. The complete SNe Monte Carlo dataset 

Monte Carlo supernova "data" points were generated 
using the red shift distributions and magnitude uncer- 
tainties specified above. For each supernova from each 
simulated experiment, the appropriate red shift distribu- 
tion was used for the random generation of a red shift 
z. (To be specific, the red shift probability distribution 
Pi{z) for each experiment i was first normalized so that 
its maximum value was 1. Then, two random numbers, 
r\ and T2, were chosen from a uniform distribution be- 
tween and 1. If the probability Pi{z max ri) at red shift 
z max r i was greater than r2, then z mal r 1 was taken to 
be the red shift of the supernova. Otherwise, the two 
numbers were discarded and the process repeated. Here, 
Zmax, the maximum allowed supernova red shift, was set 
to 2.0 for each experiment.) Next, the apparent magni- 
tude m^ ld ^ (z) expected for a supernova at red shift z was 
computed for the fiducial model used here: Q m = 0.31, 
Wo = — 1, Wi = 0. The magnitude uncertainty Um^\z) 
for that supernova was computed as described in the pre- 
vious subsection. Finally, the simulated apparent mag- 
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therefore taken to be 
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FIG. 5: The simulated SN la dataset. The points with error- 
bars are the simulated supernovae, while the solid line is the 
fiducial model. 



r TT,(theory) . r TT,(eff) 



(13) 



Note that (|13fl is a more conservative expression for the 
variance than that given by [8^ , since there is an "extra" 
factor of f s ky in the denominator of Var(TT). The in- 
formation "lost" by dividing Var(TT) by f s k y is actually 
stored in the off-diagonal components of the covariance 
matrix (which we ignored), because nearby C;s are anti- 
correlated |70j . However, we found that this extra f s k y 
had only a negligible effect on our forecast constraints. 

We used a similar form for the uncertainty in the E- 
mode polarization power spectrum Cf E . Since WMAP 
did not include the C EE s in the likelihood analysis of their 
first year data, effective corrections to the noise and sky 
fraction were not computed. We used f s k y = 0.85 and 



AC 



EE,(eff) 



EE 



nitude rrS MC ^ was computed by adding to m^ ld \z) a 
random number chosen from a Gaussian distribution, 
with mean zero and standard deviation o~m ^ ■ These 
Monte Carlo simulated magnitudes mS MC \ along with 
their corresponding red shifts and magnitude uncertain- 
ties, formed the simulated dataset used here. 



Table IIVI lists the numbers and red shift distribu- 
tions of the supernovae simulated, and Fig. plots 
the magnitude-red shift relation of the simulated "data" . 
Note that, since the major supernova surveys planned 
over the next several years are ground-based, the simu- 
lated dataset contains very few supernovae with red shift 
greater than about 1.2. 



CMB Simulation 



CMB covariance matrix 



In the following discussion, we use the convention of 
[89| that Ci = 1(1 + l)C//(27r) for any power spectrum 
C/. We computed variances in the C;s using the expres- 
sions given in Jjjj and [T^ for the diagonal components of 
the covariance matrices. The noise parameters A/" ; TT and 
Mp E were defined to be 1/8 of those used by WMAP, 
in order to simulate the effects of using eight years of 
WMAP data. WMAP noise parameters were taken from 
their data tables, which were described in [7l| and 
and provided along with the WMAP likelihood code de- 
scribed in The effective noise parameter Af, 
as a fraction of the original noise parameter J\fF T , and 

the effective sky fraction f^J y (I), were defined as in [Tof . 
The final expression for the uncertainties in the Cf T s was 



^Var(EE) 



nEE,{theory) 
'I 



+ M EE 



(14) 



Finally, uncertainties for the cross-power spectrum Cf E 



were taken from Eq. (10) in 

rTE 

' sky 



using f sky = 0.85 and 



rj-\ if £ 

f .,' 11 — 4^ as described in that reference. 



1.14 



Var(TE) 



{C TT +M TT ){C EE +Af EE ) + iC TE )2 



{21 + l)fskyf sky 



TE,eff 



(15) 

Even in the absence of a sky cut, correlations exist 
among the three power spectra Cf T , Cf E , and Cf E for 
each I |33j. As above, we keep the "extra" factor of f s k y 
in the denominators of the expressions given in that ref- 
erence. For a given multipole I, the covariance matrix 
C X Y, where X,Y = TT, TE, or EE, is given by, 



J TT,TE 



.(0 

J TT,EE 



J TE,EE 



Var(X) 

(l + 2) fsky 

(Cf E f 



G + 2) fsky 

{C EE +M EE )Cf E 

(I + \)fsky 



(16) 
(17) 

(18) 

(19) 



Large-/ correlations are also expected to exist between 
the CMB power spectra and large scale structure, due 
to the Sunyaev-Zeldovich effe ct [9(1 l9ll | and the gravita- 
tional lensing of the CMB [92ll93ll. However, these effects 
are too small to detect today |94l I05L |9(| . The analysis 
presented here neglected such correlations. In particular, 
we assumed the CMB \ 2 function to be independent of 
the red shifts z s of weak lensing sources. 
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2. Simulation 

A set of simulated TT, TE, and EE CMB power spec- 
tra was generated by means of a Monte Carlo simulation. 
For the purposes of this simulation, it was assumed that 
the WMAP project Q would be extended to eight years. 
Data from the Planck mission, as described in [93 and 
|f)8j . were not simulated for this work; the analysis of 
Planck data will be complicated by nonlinear effects such 
as the gravitational lensing of the power spectra. 

The WMAP-8 data set described here contained eight 
years of simulated WMAP data. We simulated TT data 
up to a maximum multipole of lj nax = 900, and TE and 
EE data up to l^ ax = 512. As detailed above, we divided 
by eight the tabulated WMAP noise parameters from one 
year of data. Given the CMBFAST-generated fiducial 



TE,(fid) 

l 



and C 



EE, (fid) 
I 



as well 



power spectra Q^ T ^ ld ^ £ 
as our simulated noise parameters, the power spectrum 
covariance matrices H16ll9j) were computed. For each Z, 

the inverse covariance matrix C^' was diagonalized, 
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FIG. 6: WMAP-8 simulated TT power spectrum. The points 
with errorbars are the simulated "data", while the solid line 
is the fiducial model. 



D« = U (/rl C (/rl U (i) = 



WW 







1 

(Suf) 2 




(Su 



(20) 

using the matrix XjW of eigenvectors of CJW . Next, 
U"' was used to rotate the vector of fiducial power spec- 
tra into the eigenbasis, 



C 



TT,(fid) 
I 



c 



TE,(fid) 



(21) 



c. 



EE, (fid) 
I 



The Monte Carlo simulated power spectra were 
generated in the eigenbasis, using random num- 
bers n, f2, and r3 chosen from a Gaussian 
distribution with zero mean and unit variance: 



1, (MC) 

H 

2, (MC) 

h 

3, (MC) 



T\ Su\ 

r 2 Suf 



r 3 Suf. 



(22) 
(23) 
(24) 



Finally, the Monte Carlo power spectra were rotated back 
into the standard basis, 



/ gTT,(MC) 
I r EE,(MC) 



/ u^ MC) 

t 2,(MC) 
y u 3.(MC) 



(25) 



These final simulated power spectra, qJT,(MC) ^ 

n TE,(MC) j n EE,(MG) , . 

C t v , and C ; A , were the ones used m 
our analysis. Note that this simulation method ignores 
non-Gaussianities in the distributions of observed C; val- 
ues about the underlying model; see 89] for a discussion 
of such effects. The simulated TT, TE, and EE power 
spectra are shown in Figs. ElEI and [SI respectively. 



O 




FIG. 7: WMAP-8 simulated TE power spectrum. The points 
with errorbars are the simulated "data", while the solid line 
is the fiducial model. 



D. Weak Lensing Simulation 

1. Convergence power spectrum 

In the Limber approxi mation, t he convergence power 
spectrum is given by |99l llOOl Il0l| : 



I ^ U til 



9 2 {X) 



P. 



3D 



sinjr(x) 



X)dx, (26) 



where P^d is the 31? nonlinear power spectrum of the 
matter density fluctuation, S; a(x) is the scale factor; 
and swkX = K~ x / 2 sin(X 1 / 2 x) is the comoving angu- 
lar diameter distance to \ (for the spatially flat universe 
used in this analysis, this reduces to x). The weighting 



11 



0.02 r 
0.015 - 
2= 0.01 




I 



FIG. 8: WMAP-8 simulated EE power spectrum. The points 
with errorbars are the simulated "data", while the solid line 
is the fiducial model. 
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FIG. 9: Monte Carlo simulated convergence power spectrum. 
The points with error boxes are the simulated "data" for a 

survey with f s k y = 0.7, n « 6.6 x 10 8 sr _1 and (j^t) 1 ^ 2 ~ 0.4, 
while the solid line is the fiducial model. 



function g(x) is the source-averaged distance ratio given 
by 



i(x) 



XH n(x' f* K{X '- X) d X ', (27) 
am K {x) 



where n(x(z)) is the source red shift distribution normal- 
ized by J dz n(z) = 1. We assume that all the sources 
are at a characteristic red shift z s , so that lf2*7|) reduces 
to 



sin-xiXs - x) 
shiidxs) 



(28) 



For weak lensing calcu lations, we use the standard BBKS 
trans fer fu nction |l02j | , and the analytic approximation of 
Ref. |103| for the g rowt h factor. We use the mapping pro- 
cedure HALOFIT 104J to calculate the non-linear power 
spectrum. 



2. The Weak Lensing \ 2 function 
For t he we ak lensing spectrum, the uncertainty is given 

by: mEml 



(21 + I) f sky 



P, K 



(lint 2 ) ' 



(29) 



where f s k y = O 2 7r/129600 is the fraction of the sky 
covered by a survey of dimension O in degrees, and 

1/2 

(lint) ~ 0-4 is the intrinsic ellipticity of galaxies. We 

0.7, in the same 
. We used an av- 



consider a reference survey with f ^,, = 
range as the Pan-STARRS survey (105 
erage galaxy number density of n ~ 6.6 x 10 8 sr _1 . \ 2 was 



computed between multipoles l m in and l m ax- For l m i„, 
we took the fundamental mode approximation: 



In 



360 deg 




fsky 



(30) 



i.e. we considered only lensing modes for which at least 
one wavelength can fit inside the survey area. We used 
lmax = 3000, since on smaller scales the nonlinear ap- 
proximation HALOFIT to the nonlinear power spectrum 
may not be valid. 



3. Simulation 

The procedure described above was used to generate 
fiducial convergence power spectra p^^ ld \ We used 
the cosmological parameter fiducial values of Sec. Ill Al 
We then calculated the uncertainties on the convergence 
power spectra using Eq. I|29|l . The Monte-Carlo simu- 
lated convergence power spectra are generated using 



pK,MC _ pK,(fia 



r5P? 



(31) 



where r is randomly choosen from a Gaussian distribu- 
tion of standard deviation 1 62]. The simulated conver- 
gence power spectrum for the reference survey and the 
fiducial model are shown in Fig. 



IV. RESULTS AND DISCUSSION 
A. Goals 

As we have said, the most important question in the 
study of dark energy today is whether the dark energy is 
or is not a cosmological constant. If dark energy is shown 
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FIG. 10: Current constraints on dynamical dark energy pa- 
rameterized by I^J, using first-year WMAP data, the SNe la 
"gold set" from [g, and the current SDSS data. The la, 2a, 
and 3a contours are shown. The region of parameter space ex- 
cluded by the wo + wi < prior is filled with upwards-sloping 
stripes. 
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FIG. 11: Same as Fig. 1101 except that the dark energy param- 
eterization © has been used. The region of parameter space 
excluded by the wo + w a < prior is filled with upwards- 
sloping stripes. 



not to be a cosmological constant, the next questions 
that arise are again qualitative: Is the equation of state 
greater or less than —1? Does it vary with time? The 
actual values of u>(0) and w'(Q) will become important in 
the more distant future, when theoretical models for the 
form of w(z) are available. We look for constraints on 
the equation of state parameters, not because we believe 
that any one of the parameterizations Q [21 E El is a 
theoretically favored model for w(z), but because we wish 
to address the above qualitative questions about dark 
energy. Thus, we emphasize that our goal is to answer 
these qualitative questions in a careful, parameterization- 
independent way, rather than to find specific values for 
w(0) and w'(0) within a particular parameterization. 



B. Analysis of Current Data 

Keeping these goals in mind, we used \ 2 minimization 
to analyze current data from the CMB power spectrum 
, the SN la "gold set" , and the galaxy power spec- 
trum 12]. We began by using to parameterize the 
dark energy equation of state. The contours obtained 
are plotted in Fig. with grid spacings A Wo = 0.06 
and A Wl — 0.15. The resulting la and 2a constraints 
on the dark energy parameters are wq = — 1.381q 20^0 55 
and wi = 1.2±g;fg±5;gf , with a best-fit \ 2 value of 1611. 
We consider our 2a contours to be more reliable than our 
la contours, since the latter occupy an area in Fig. 1101 
only a few times the area of a grid square. 

ACDM models are not conclusively ruled out by our 
results. Compared to the best-fit model in Fig. ^| the 
ACDM model has a x 2 value that is higher by 5. How- 
ever, the cosmological constant is a "simpler" model of 
dark energy than the w(z) parameterizations considered 
here, in the sense that the ACDM model has two fewer 



variable parameters than JTJ [3 EJ. Thus the ACDM 
model has x 2 /d-o.f. = 1616/1514 = 1.0676, while the 
best-fit model from Fig. [TOlhas x 2 /d.o.f. = 1611/1512 = 
1.0656. In the approximation that uncertainties in the 
data points are Gaussian, these correspond to probabil- 
ities of Pkcdm = 0.0336 and P best ~ftt = 0.0378. For 
comparison, a point on the edge of the 2a contour in Fig. 
Unihas x 2 /d-o.f. = 1615/1512 = 1.0682 and probability 
0.0324. In this sense, the ACDM model is more proba- 
ble than some points within the 2a contour. Although 
the data favor u> < — 1 and wi > 0, the cosmological 
constant is not decisively ruled out as a model of dark 
energy. 

Moreover, we have not considered systematic uncer- 
tainties in the supernova data, which can degrade con- 
straints on the equation of state. As an optimistic es- 
timate of the effects that such a systematic will have 
on dark energy constraints, we assumed an uncorrelated 
systematic uncertainty 5m = 0.04, in red shift bins 
Az = 0.1, and reanalyzed the data. The added uncer- 
tainty caused the 2a contour to broaden to include the 
cosmological constant, while the best-fit model remained 
the same. 

It is troubling that the best-fit model in Fig. ^|is so 
close to the boundary wq + w\ = of the dark energy 
parameter space. This is despite the fact that |0J, out 
of the three 2-parameter equations of state illustrated 
in Fig. ^ covers the greatest range of parameters in 
the (z, w(z)) plane. Parametrization J2J, widely used 
in the literature, is more restrictive than |0J, and there- 
fore should have more problems with this boundary. We 
repeated our analysis using J3J, obtaining the \ 2 con- 
tours shown in Fig. II II with grid spacings = 0.1 and 
A Wa = 0.25. The resulting 2a bounds are w = — l-3l .34 
and w a — 1.25^217, with a best-fit \ 2 value of 1613. 
Thus we see that the allowed region shifts significantly 
when we switch to parameterization J2Jl, with u/(0) < 
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FIG. 12: Comparison of the two dark energy parameteri- 
zations Contours corresponding to |QJ| are filled in 

with vertical stripes (x 2 = 1615 contour) or upwards-sloping 
stripes (x 2 = 1620 contour), and the best-fit model is marked 
by an "X". Contours corresponding to @ are filled in with 
horizontal stripes (x 2 = 1615 contour) or downwards-sloping 
stripes (x 2 = 1620 contour), and the best-fit model is marked 
by an "O". Note that the x 2 = 1615 and x 2 = 1620 contours 
are the 2a and 3<r contours, respectively, of parameterization 
@l . The thick and thin dashed lines correspond to the lines 
wo + wi = and Wo + w a = 0, respectively. 



models falling within the 2a contours. Also, the mini- 
mum x 2 value goes up by 1.8. 

These shifts are better understood when the \ 2 func- 
tions corresponding to the two parameterizations are 
plotted together in the (wq, w') plane (with w 1 = w(l) — 
w(0)), as in Fig. Note, first of all, that the two sets 
of x 2 contours would be nearly identical if the region of 
parameter space between the dotted lines were removed. 
This region corresponds to the slice of the (z, w(z)) plane 
between the dashed and solid lines in Fig. ^ which is al- 
lowed by parameterization (0J but not by Also note 
that the best-fitting point (wq — —1.38, w' = 1.2) of 
parameterization Q), marked by an "X" on the plot, is 
outside the range allowed by parameterization (0) . Since 
the minimum x 2 found using J2J) is higher, the 2a con- 
tour is at a higher x 2 1 meaning that it now encloses some 
points with w a < 0. Thus, the difference between the 
constraints found using (0J and (J2J) is directly related to 
the parameterization, and particularly, to the high red 
shift w(z) constraint discussed in Sec. Iff Al 

To summarize, switching parameterizations from (0J 
to (J2J) moves the boundary downwards, pushing the con- 
tours in the direction of decreasing w' . Since the contours 
in Fig. HOI are still near the boundary, a parameterization 
that allows larger values of w' may allow the contours to 
shift even more in that direction. It is possible that such 
a parameterization would favor wq < — 1 and w' > 
to an even larger degree, leaving the ACDM model less 
favored. 

Comparison with the published literature shows that 
our results are consistent with others obtained using var- 
ious combinations of cosmological probes. Table IVl lists 



several recent analyses, along with the parameters and 
w(z) parameterizations used (parameters w' , w a , and wi 
imply parameterizations (Q, J2J, and QJ, respectively), 
the priors, and the resulting w(z) constraints. When 
w(z) constraints using the same equation of state pa- 
rameterizations are compared, our results are consistent 
with those in Table IVl at the 2a level. In addition, the re- 
sults of |8l ll07l |. obtained using parameterization and 
SN la data, are consistent with our constraints using ei- 
ther of the parameterizations EJ ■ It is reassuring that 
our x 2 minimization procedure, which handles parame- 
ter degeneracies differently than marginalization, obtains 
consistent results. 

Several papers list constraints on w(z) at specific val- 
ues of shown in Table IVT1 It is evident from the ta- 
ble that u>(0.3) is well-constrained and parameterization- 
independent. Our results for w(0.3) are in excellent 
agreement with those of |56| . even when their three- 
parameter equation of state (J3J) is used. Meanwhile, the 
equation of state at z = 1 is less well constrained, though 
once again our results agree with those of 56J at the 2a 
level. Parametrization Q evidently prefers higher val- 
ues of w(l) than does One would expect the addi- 
tion of the third parameter in © to bridge the gap be- 
tween constraints using @ and I0J. However, our best-fit 
models using (@J have w(0) « —1.4, u>(0.3) « —1, and 
w(l) w —0.2. Even © cannot reproduce all three of 
these features if one also imposes the CMB constraint 
w(z) < at z > 1, which implies wq + w a + < 0. 
That is, setting w(0) = -1.38 and w(l) = -0.18 from 
our best-fit model, and imposing the requirement that 
lim z _ >00 w(z) < 0, implies that w(0.3) > —0.70 in pa- 
rameterization Q , which is inconsistent with constraints 
on w(0.3) discussed above. Since our best w models are 
qualitatively different from those allowed using the pa- 
rameterizations © and ©, there is still no discrepancy 
between our results and those of (56J for w(l.0) Finally, 
recall that the two-parameter equations of state Q 
0} relate w(z) and its derivative, at low red shifts, to 
the equation of state in the large red shift limit (see Ta- 
ble [J. Our constraints on the large red shift value of 
the equ ation of state, lirnz^oo w(z), agree closely with 
those of |110| . These results appear to be independent of 
parameterization; analyses with three different parame- 
terizations agree that w(z) at large red shifts is slightly 
less than zero. 

Our uncertainties in the dark energy parameters are 
mostly in agreement with similar analyses that include 
CMB and SN la data. Our 2a uncertainties in wq and 
WQ + w a = lim^oo w(z) are appr oximately twice as large 
as the 68.3% constraints of |110| . just as expected. Also, 
our 2a constraints on wq ar e about twice as large as 
the la uncertainties of |l06j . Their looser bounds on 
w'(0) can be attributed to the fact that they use a four- 
parameter equation of state, and marginalize over param- 
eters other than those reported. Meanwhile, comparison 
of our results with |5fj| provides an example of parame- 
terization effects on the dark energy constraints. Their 
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Ref. 


probes included 


parameters varied 


priors 


dark energy constraints 


bl 


SN la 


wo, w', Q m 


Q m = 0.27 ±0.04 


w = -l.31t_.__, 

' 1 zio+O-81 

w = 1.48I 090 




rtMfD CAT T~ TD/U\ 

CMn, oiv la, r\K), 
bias, Ly-a 


"vanilla" , wo , w a 


— ^ Cl o 
T < U.O 


nQ81 +IJ.W3W.3«4+0.l.8« 
WO — U.y<51_ 193 _o. 373-0. 521 > 
,,, _ n nr+0.65+1.13+1.38 
w a — u.uo_ g3 _ 1 g2_2. 88 




kjLvld, ji\ la, i \ j 


"vanilla" , 

4 w(z) parameters 


h 79-1-0 08 O^O 1 

0.014 < io b < 0.04, < t < 1, 
6 < n < ~\ A 

KJ • \J ^ It's -L • 


Wo — — 1.4d_ 38 , 

dw I 1 ri + 1.0 

dz \ z=0 ~ 1 - U -0.8 


107] 


SN la 


w , w', Q m 


d%c } uU' 2 /c = 1.716 ± 0.062, 
d ,\ nD \ , -0.51 ±0.11 

din a | * = n.1K 


-1.05 < w < -0.29, 
-1.89 < w < 0.05 a 


1071 
I 1 " "J 


SN la 


i ,in n\ O 

wo, w , Jim 


o — n 97 4- n oa 


— 1 3Q <r j/in <r —0 95 
-2.61 < w' < 1.49 a 


i no 1 

108] 


CAT 7 „ 

oiv la 


Q m , 3 w(z) parameters 


.(dec) 1/2 / 1 71(1 J-n 157 

a AC u>m jc — I. ( 1U ± U.lof, 
u b = 0.024, u m = 0.14 ± 0.02 


,T<\ 1 OQ7+0-016 
u 'lo<*<0.414 — -0.056' 
, 7 .| __ 22 o+0.070 
W \0A14<z<l — U ' ZZ9 -0.117' 
,7,1 — f) 149+0-05 1 
W \l<z<1.755 — U ' i ^ z -0.033 


[ioi 


SN la 


w , w a 


a m = 0.27 


WO = -1.38 ±0.21, 
w a = 2.78 ± 1.32 


109] 


SN la 


w , w a 


VL m = 0.31 


w = -1.35 ±0.24, 
Wa = 2.32 ± 1.56 


i — 

109] 


SN la 


WO, Wa 


fl m = 0.35 


wo = -1.3 ±0.29, 
w a = 1.64 ± 1.95 


110] 


CMB, SN la, X-ray 


"vanilla", bias, wo, w a 


bias = 0.824 ± 0.089, 
L0 b = 0.0214 ±0.0020, 
h = 0.72 ±0.08 


wo = -1.16__;_|, 

wo ± w a = -0.05__;_?. 


lud 


CMB, SN la, X-ray 


"vanilla", bias, Qk, wo, w a 


bias = 0.824 ± 0.089, 
uj b = 0.0214 ± 0.0020, 
h = 0.72 ±0.08 


w = -1.14t_;__, 

wo ± w a = -0.09_°;__ 


110] 


CMB, SN la, X-ray 


"vanilla", bias, 
3 w(z) parameters 


bias = 0.824 ± 0.089, 
LOb = 0.0214 ± 0.0020, 
h = 0.72 ±0.08 


wo = -1.231S:S, 
lim z ^ 00 w(2) = -0.12±g;:£ 


AU, 

MI 

PJS 


CMB, SN la, P(k) 


"vanilla", wo, w± 


for priors see Table ITTI 


1 QS+H.3U 

w - -1.38_ 06S , 

1 on+0-64 
wi = l.20_ 106 

(2a constraints) 


AU, 

MT 

PJS 


CMB, P(k), SN la 


"vanilla", wo, wi 


for priors see Table ITTI 


w = -1.36t^, 

i on+0.60 
Wi = 1.20I 188 

(2a constraints) 


AU, 
MI, 
PJS 


CMB, SN la, P(k) 


"vanilla", wo, w a 


for priors see Table ITTI 


„, _ i o+U.39 

w a = 1.25t°;f 7 
(2a constraints) 



"Bounds taken from plot of 68% probability contour 

TABLE V: Dynamical dark energy constraints from current data. The "vanilla" parameter space is spanned by the six param- 
eters (h, oj m , w;,, t, A, n s ). The cosmological probes used are: bias = SDSS galaxy bias, CMB = WMAP CMB power spectra, 
Ly-a — SDSS Lyman-a, P(k) = SDSS galaxy power spectrum, SN la = SN la "gold set", X-ray — Chandra X-ray clusters. 



95.5% upper and lower bounds on w are nearly the same 
as our 2(7 bounds, and their 95.5% lower bound on w a is 
somewhat tighter than ours. However, our upper bound 
on w a is smaller than theirs by a factor of three, due to 
the fact that our best-fit model is much nearer to the 
boundary wq + w a — of parameter space. Still, the 
boundary has some effect on their 95.5%, and 99.86% 
probability bounds on w a . Although their 68.32% upper 
and lower bounds are nearly the same, their 99.86% up- 
per bound is smaller than the lower bound by a factor 
of two. Thus, it would be interesting to see their analy- 
sis repeated with parameterizaton (J3J|, which has a less 



restrictive boundary. 

The effects of systematic uncertainties can be esti- 
mated by comparing the analysis procedures used by the 
studies listed Table Ivl as well as the locations of the re- 
sulting contours in parameter space. Ref. finds 68% 
an d 95% contours centered near the ACDM model, while 
d HoE Eil HHl [nil as well as our analysis find 
contours that lie mostly in the wo < —1, w'(0) > re- 
gion. Even when the analysis of |56J is repeated without 
the Lyman-a and galaxy bias data sets, th eir p referred 
models are still close to the ACDM model |ll2j . Since 
both studies use the standard SN la likelihood, the differ- 
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z 


Ref. 


Parameterization 


w(z^) constraint 


0.3 


loo] 


i w% 


-1.011I .2i5 [9o.57oJ 


150] 




loll 


— u.yoi_ 249 (,yo.0/oj 


ATT T\ TT TJ> TO 

AU, Ml, rja 


H7T 


-i-OS-o.as ( 2<J ) 


ATT Tl |T T~» TO 

AU, Ml, PJb 


m 


- 1 - 01 -0.29 ( 2(J ) 


1 


156] 


mi 


-i-OO-o.es (95.5%) 


156] 


iF?l 

ll.il 


— l.Uo_ 58 (yo.o/oj 


ATT Tl TT T> TO 

AU , Ml, PJb 


;m 


-0.18-o.8o ( 2ct ) 


ATT T\TT T3 TO 

AU , Ml, PJo 


F7I 


-0.68_ 74 (2cr) 


oo 


rno] 


© 


-0.05+H 9 7 


[110] 


3 w(z) params. 


-O.i2± u ;^ 


AU, MI, PJS 




-0.18+-- (2a) 


AU, MI, PJS 




-0.05--- (2a) 



TABLE VI: Constraints on w(z) at several redshifts. 



Analysis 


Par am. 


wo (to 2a) 


Wi or 7j; a (to 2cr) 


standard 


a 


-1.0»_o.55 


"1=1.21^ 


standard 




i O+0.39 
1 -°-0.34 


w a = L25±£# 


CMB I > 20 




-1 63 +u ' 43 

J -°°-0.65 


-i = i-5ir 6 , 


wo + 2wi > 0, 
CMB I > 20 




1 9 t-+U.40 


«ll = O.6I1 4 , 



TABLE VII: Consistency checks of our current w(z) con- 
straints. All constraints shown are at the 2<r level. The "stan- 
dard" analyses above use WMAP, SDSS, and SN la gold set 
data, with priors as shown in Table ITT1 



ence between their result and o urs m ust come from cither 
their CMB likelihood f uncti on |l!3j or th eir galaxy power 
spectrum analysis [HE Ill4j . Ref. claims that the 

WMAP likelihood approximation is inaccurate at low I, 
and that the WMAP foreground removal procedure may 
lead to a suppression of low-/ power. Their likelihood 
function is designed to correct these problems. Mean- 
while, their galaxy power spectrum analysis uses SDSS 
measurements up to k = 0.2/i/Mpc, a range which ex- 
tends into the nonlinear regime. The combined effect 
of these two changes is to move the contours very close 
to the ACDM model, which is disfavored by A\ 2 = 3 in 
our analysis with the same w(z) parameterization. Thus, 
differences in the CMB and galaxy power spectrum like- 
lihood functions lead to a shift of over la in the x 2 con- 
tours. 

We attempted to reproduce this shift by using only 
the multipoles / > 20 in our CMB analysis. If the low 
multipolcs were responsible for the shift, then neglecting 
them would enlarge the \ 2 contours and move them to- 
wards the ACDM model. However, as shown in the third 
column of Table IVIII the allowed region actually moves 
away from that of 56], although it does broaden some- 
what. Not only does this test fail to explain the difference 
between the regions of parameter space preferred by j5|| 
and our results, it also demonstrates that the CMB anal- 
ysis is sensitive to the low / region of the power spectrum. 
Since there is an ongoing debate about the proper han- 



dling of the low multipoles 113], this sensitivity to that 
region of the power spectrum is worrisome. 

Recall that switching dark energy parameterizations 
can change the uncertainties in w(0) and u/(0). In order 
to test the effects of this parameterization dependence in 
combination with the low / sensitivity, we imposed the 
prior constraint wa + 2wi = wq + 2w' > on the anal- 
ysis with CMB I > 20. This allowed us to estimate the 
constraint on w a that we would have obtained if we had 
repeated the analysis with parameterization @ ■ The re- 
sulting 2a constraints, shown in the fourth row of Table 
IVIII are very similar to those found by the standard anal- 
ysis using parameterization |2"|l. They have not moved 
appreciably towards the ACDM model, though the uncer- 
tainties have grown. The slight shift away from ACDM 
seen without the prior wq + 2w' > has been hidden by 
the imposition of that prior. This implies that the pa- 
rameterization dependence of the results conceals some 
of their sensitivity to low I CMB data. Once again, we 
are unable to reproduce the shift in the location of the 
equation of state contours. The explanation may have 
to do with the details of the CMB likelihood function 
used by [5||, or with their galaxy power spectrum anal- 
ysis. These differences between the likelihood functions 
should be examined critically. 

Similarly, the SN la likelihood function used by |107| 
causes the dark energy constraints to shift relative to 
those in 8] . The standard supernova analysis Q predicts 
-1.59 < w < -1-09 and 0.58 < w' < 2.29.' Using the 
same dark energy para meterization and the same prior 
fl m = 0.27±0.04, QJ22I obtains -1.39 < w Q < -0.25 and 
—2.61 < w' < 1.49 . The |107| analysis uses supernova 
flux averaging |l!5l | to handle s ystem atic effects due to 
the weak lensing of supernovae |ll6j . The flux averag- 
ing method assumes that uncertainties in SN la fluxes, 
rather than in magnitudes, are Gaussian. If uncertain- 
ties ar e act ually Gaussian in magnitude, this introduces 
a bias |115| : conversely, if uncertainties are Gaussian in 
flux, then the standard SN la analysis is biased. The 
net effect of the flux averaged likelihood function is to 
weaken constraints on wo and w' , and to shift them both 
by about la towards the ACDM model. 

Another SN la systematic effect, not included in any 
of the analyses in Table [V] is the dimming of supernovae 
by astrophysical dust. Ref. |117| considers several types 
of intergalactic dust, and constrains supernova dust dim- 
ming to be less than 0.2 magnitudes. (For an extensive 
discussion of SN la systematics and their effects on the 
cosmological parameters, see |76|.) The upper bound on 
dust dimming is not much less than the mean SN la mag- 
nitude uncertainty in the current data set, which suggests 
that dust dimming can cause a significant bias in dark 
energy constraints. Since the dimming of SNe la pro- 
vides evidence for accelerating cosmological expansion, 
additional dimming due to dust will appear to exagger- 
ate this acceleration. Ignoring systematic uncertanties 
due to dust dimming will cause an analysis to underes- 
timate w, since more negative values of the equation of 
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state lead to a more rapid acceleration. 

To summarize, current data are not precise enough to 
address whether or not the dark energy is a cosmological 
constant. A careful study of the current data, and com- 
parisons with other recent analyses, reveal several obsta- 
cles to a satisfactory understanding of the dark energy. 
The degeneracy between wq and w\ remains unresolved, 
and the allowed 2a ranges for wq and w\ are at least 
-1.93 < w < -1-08 and 0.14 < w x < 1.84. Further- 
more, this result is dependent on the parameterization 
chosen for the dark energy equation of state. In particu- 
lar, problems arise when the la and 2a contours lie near 
a boundary of the dark energy parameter space, as they 
do at present. More generally, uncertainties in the equa- 
tion of state parameters depends on the location of the 
best-fit model in the parameter space. Aside from the pa- 
rameterization dependence issues, analyses that use the 
same data, but different treatments of the systematics, 
find best-fit models differing by ~ la. The combination 
of parameterization dependence and systematic effects is 
enough to shift the best-fit model by more than 2a. Thus 
we cannot rule out the ACDM model at present. 



C. Forecasts using Simulated Data 

Next we ask, when data from future probes are avail- 
able, how close we can come to achieving our goal of 
distinguising between a cosmological constant and other 
forms of dark energy. The dark energy parameter space 
is continuous, with nonconstant w(z) models arbitrarily 
close to the ACDM model. Unless we are lucky and find 
a dark energy significantly different from the cosmologi- 
cal constant, the best we can do is to distinguish between 
models separated by some minimum distance in param- 
eter space. 

Thus, we begin by identifying a set of reasonable mile- 
stones against which progress in w(z) constraints may 
be measured. Since we do not know ahead of time which 
best-fit model will be found, we will say that an exper- 
iment can distinguish between two points A and B in 
parameter space if, regardless of the location of the best- 
fit model, the 2a contour around that best-fit model ex- 
cludes either A or B (or both). 

Our first milestone is to distinguish between a cosmo- 
logical constant and a "dark energy" whose equation of 
state at z = 1 is w(l) = 0. This goal is motivated by 
current observations, rather than by theoretical consid- 
erations; several of the analyses listed in Table in- 
cluding our own, favor w ~ at red shifts of order unity. 
We can distinguish between these two types of dark en- 
ergy by ruling out cither w(l) = —1 (ACDM model) or 
w(l) = 0. In the worst case, the best-fit equation of state 
will have w(l) = —0.5, meaning that the 2a uncertainty 
in w(l) must be brought below 0.5 in order to rule out 
one of the two types of dark energy. Thus, the first mile- 
stone is decrease the uncertainty in w at z = 1 to less 
than 0.5. 



As our second milestone, we would like to distinguish 
between the ACD M mo del and tracker quintessence |l!8j 
or tracker SUGRA [H3 models with w > -0.8, wi > 0. 
These are interesting, partly because they represent a 
large class of well established, theoretically motivated 
models, and also because they make definite predictions 
that are clearly distinct from the ACDM model in pa- 
rameter space. It is not clear whether such models are 
ruled out by current data; our analyses exclude them to 
greater than 2a, but they are within the 95% probability 
contour of |56l |. In the future, we can compare an ar- 
bitrary best-fit model to ACDM and tracker models by 
comparing their values of w(z*), where is the red shift 
at which the uncertainty in w(z) is a minimum. The 
ACDM model predicts w(z Sf ) = — 1, while trackers have 
w(z*) > —0.8. In the worst case, the future best-fit model 
will have w(z*) ~ —0.9. Therefore, in order to reach the 
second milestone, we must reduce the 2a uncertainty in 
w(z*) to about 0.1. 

Finally, our third milestone is to distinguish between 
ACDM and quintessence models that are near to it. Wc 
do not have in mind any particular class of theoretical 
models, so "near" is not well-defined. As a reasonable 
third milestone, we consider aiming to distinguish ACDM 
from models with either \wq + 1| 0.05 and w\ — 0, or 
wo = —1 and |iui| > 0.05. For this we will need both 
Aw (2a) < 0.025 and A Wl (2a) < 0.025, that is, a Wo ~ 
Cu>i ~ 0.01. If the ACDM model were to be the preferred 
model even when the equation of state uncertainties were 
reduced to this level, then it would probably be time to 
abandon our hopes of using w(z) to study the nature of 
dark energy. 

Now that Sec. IIV Bl has pointed out several pitfalls 
in the analyses of current data, we make the resonably 
optimistic assumption that the likelihood functions asso- 
ciated with the CMB, SNe la, and cosmic shear will be 
well understood by the end of the decade. We assume 
that these new, accurate likelihood functions will give no 
less cosmological information than the likelihood func- 
tions discussed in Sections llff BlIjTTCl and llll Dl This as- 
sumption is optimisti c; the conservative CMB foreground 
removal procedure of |l!3j | , as well as the inclusion of SN 
la dust dimming effects, should increase uncertainties in 
the dark energy parameters. Moreover, for weak lens- 
ing, several systematic effects (e.g. intrinsic alignments 
of source galaxies, selection biases and residuals from the 
PSF correction) need to be better understood and tightly 
control led i n order for this probe to achieve its full po- 
tential [l2Cj . 

Also, our forecasts are based on a ACDM fiducial 
model, since (wq = —1, w'(0) = 0) is far from both the 
parameter space boundaries wq+wi = and wo+w a = 0. 
From now on, we limit our study to the parameterization 
alone. Assuming that our x 2 contours stay away from 
the boundaries, we can compare constraints between pa- 
rameterizations (0J and J2J using the "rule of thumb" 
&wi = &w' ~ cr iUa/2, where the average low red shift 
derivative w' = w(l) — w(0). The factor of 1/2 means 
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that, even when parameter space boundaries are unim- 
portant, the uncertainty in the derivative u/(0) is param- 
eterization dependent. 

Dark energy constraints provided by eight years of sim- 
ulated WMAP data alone are shown in Fig. with grid 
spacings A WQ — 0.3 and A Wl — 0.6. Models in the upper 
right-hand section of Fig. 1131 are ruled out by the prior 
constraint w + wi < 0. Meanwhile, the requirement that 
h < 1.1 eliminates models in the lower left-hand corner 
of the plot. Therefore, to 2a, virtually no constraints 
are imposed on w and Wi within the ranges shown in 
Fig. Such weak CMB constraints on dynamical dark 
energy are to be expected, given the angular diameter 
distance degeneracy discussed in Sec. Ill ±31 

Given the fairly weak priors used here, the supernovae 
alone were also unable to provide interesting constraints 
on dynamical dark energy. Analysis of the supernova 
dataset shown in Fig. [5]gave deceptively tight constraints 
on wi; see Fig. ^] with grid spacings A,„ = 0.02 
and A Wl — 0.1. However, multiple Monte Carlo sim- 
ulations of SN la datasets yielded 2a constraints that 
varied widely from simulation to simulation. Three out 
of our five simulations had 2a contours extending into 
the W\ < —8 range. 

The problem, as pointed out in |l2l| . is that the dark 
energy density pde(z) exhibits qualitatively very different 
behaviors for w\ > and wi < 0. When w\ > 0, Pde{z) 
can remain a nontrivial fraction of the total energy den- 
sity of the universe up to red shifts of order unity. When 
wi < 0, pde(z) drops quickly with increasing red shift, so 
that the dark energy is important only in the very recent 
past. Distinguishing between different dark energy mod- 
els is difficult when pde(z) is very small. This explains 
the fact that several of the simulated 2a contours extend 
deep into the wi < region even though they are all 
simulations of a model with wi = 0. 

Adopting strong priors on f2 m significantly reduces 
these non-Gaussianities, as shown in Fig. 1151' top'). When 



we fix f2 m = 0.31 and repeat the analysis of five simu- 
lated data sets, we find mean 2a uncertainties on wq and 
Wx of Aw {2a) = 0.522 and Awi(2a) = 1.63. The stan- 
dard deviations in Awo(2a) and Aw\{2a) are 0.036 and 
0.21, respectively. Thus when Q m is fixed, the uncertain- 
ties vary by only a few percent from one simulation to 
another. 

However, |122| points out the perils of assuming strong 
priors in the SN la analysis. In particular, if the analysis 
of another cosmological probe is based on the assumption 
that w(z) — —1, then the constraints from that analysis 
should not be imposed as priors in a study of dynamical 
dark energy. Figs. llSf mid - ) and llSf bot - ) show the effects 
of weakening the prior constraint on fl rn to fl rn = 0.31 ± 
0.05 and fi m — 0.31 ± 0.10, respectively. Even with the 
relatively strong priors used in Fig. I15f mid). the non- 
Gaussianities have returned, and constraints on the dark 
energy parameters are weakened. Thus it is imperative 
that the supernovae be combined with another data set 
in order to provide constraints on dynamical dark energy. 

Combination of the WMAP-8 and SN la datasets im- 
proved dark energy constraints considerably, as shown in 
Fig. El w hh grid spacings A Wo = 0.08 and A Wl — 0.22. 
As with the supernova analysis, five simulated datasets 
were analyzed separately in order to determine the dark 
energy constraints. We obtained the 2a constraints 
Aw (2a) = 0.52 and Awi(2a) = 1.65. The standard de- 
viations in the 2a wq and w\ uncertainties were 0.027 and 
0.25, respectively. This implies that our forecast uncer- 
tainties are each correct to about 10%. We checked these 
constraints using a Fisher matrix calculation and found 
the 2a constraints 0.44 and 1.46 in wq and w\, respec- 
tively, consistent with our minimization results. Mean- 
while, the minimum 2a uncertainty in w(z) was found to 
be Aw(z»)(2a) = 0.19 at = 0.27. 

Our forecast uncertainty, in each of the equation of 
state parameters, is greater by a factor of 2-3 than that 
of 40] with SNe la and a prior on f2 m , as shown in Table 
IVIIII Note that our la uncertainties have been computed 
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contours are shown. The location of the fiducial model is 
Wq marked by a "+" sign. 



FIG. 15: SN la contours for the simulated dataset shown 
in Fig. assuming the prior constraints Q m = 0.31 (top), 
n m = 0.31±0.05 (mid), and tt m = 0.31±0.10 (bot). The la, 
2a, and 3a contours are shown. The location of the fiducial 
model is marked by a "+" sign. 



by dividing our 2a uncertainties by two. 

Next, the analysis was repeated with two different fidu- 
cial models, chosen to lie approximately along the line of 
degeneracy in the (wq,wi) plane. As shown in Fig. El 
the model (wq — — 1.3, w\ — 0.5) had 2a uncertainties of 
0.73 and 2.58 on wq and w\, respectively, which are sig- 
nificantly larger than those reported above. Meanwhile 
the model (wq = —0.7, w% — —0.3) had uncertainties of 



0.48 and 1.58. These are consistent with our constraints 
obtained using the ACDM fiducial model. 

Returning to our ACDM fiducial model, we repeated 
the analysis with different assumptions about the data in 
order to test the robustness of our forecast constraints. 
Table IIXI lists the new constraints obtained. The stan- 
dard deviations in our original constraints are about 10%, 
so we did not consider a modification to be significant un- 
less it changed constraints by at least 20% — 30%. First 
we checked to what extent a CMB simulation without the 
"extra" factor of f s k y in the denominator of the covari- 
ance matrix (TTTTT IT4l ITTn [TTjl ITTl m?l ITTTT ) would improve 
our constraints. The difference turns out to be negligible. 
In fact, the final constraints are relatively independent of 
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Ref. 


surveys 


parameters varied; 


forecast 




included 


priors 


constraints 


[38] 


Planck + WL 
(1000 deg 2 ) 


"vanilla", u>o, T/S; 
w' = 


a WQ 


= 0.15 


[39] 


WL + COBE 
± photo z 
+ Planck 


u> m , uJb, n s , m v , wo; 
w' = 0, 
ffW m = -064, 
flnwj = .035, 
a n = .04, 
Cm„ = -58 


a-WQ 


= 0.19 


[40] 


SNAP SNe la 


Wo, w a , 0, m ; 


Vivo 


= 0.09, 




(syst. Sm — 


<7Q,„ = 0.03, 


a & 


= 0.31 




0.0074(1 + z)) 


Q K = 






m 


SNAP SNe la 


wo, Wa, fi m ; 


<Jw 


= 0.09, 




(syst. Sm — 


Planck priors, 


&w' 


= 0.19 




0.0074(1 +z)) 


Q.K = 






m 


SNAP SNe la 


WO, W a , 


(Two 


= 0.05, 




(syst. Sm — 


Or = 


(Tw> 


= 0.11 




0.0074(1 + 2)) 
+ WL 










(f sky = 0.025) 








[41] 


Planck + WL 


"vanilla" , Wo , w a , 


(Tw 


= 0.056, 




(fsky = 1, 

tomography) 


UJ V , OLs, yHe 




= 0.087 


[41] 


Planck + WL 


"vanilla" , wq , w a , 




= 0.076, 




(fsky = 0.5, 
tomography) 


ujv, a s , yne 




= 0.11 


[41] 


4 yr. WMAP + 


"vanilla" , wo , w a , 




= 0.064, 




WL (f sky = 1, 
tomography) 


Uv, CKs, yHe 




= 0.11 


[42] 


Planck + 


"vanilla" , wo , w a , 


Cw n 


= 0.05, 




SNAP SNe la + 


(Jlv, CXs, yHe 




= 0.1 




WL (f aky = 0.5, 
tomography) 








[43] 


Planck + LSST 


LUm, ^de, f8, 


(Tw 


= 0.036, 




cluster counts 


n s , w , w a ; 




= 0.093 




and power 


linear biasing 








spectrum 
(200,000 clusters) 








[123] 


1280 SNe " 


wo, w a , fi m ; 


(Two 


= 0.27, 






Planck priors, 


(Tw> 


= 0.57 






<7Q m = 0.03, 










Ok =0 






AU, 


8 yr. WMAP + 


"vanilla", wo, Wi; 




= 0.26, 


MI, 


~ 2000 SNe 


for priors see 


(Tw\ 


= 0.82 


PJS 


(syst. Sm = 0.04) 


Table El 






AU, 


8 yr. WMAP + 


"vanilla", wo, wi; 


0~w a 


= 0.10, 


MI, 


~ 2000 SNe 


for priors see 


(Tw\ 


= 0.18 


PJS 


(syst. Sm = 0.04) 
+ WL 


Table El 







'irreducible magnitude systematic Sm' rr = 0.01 + O.O62 for 2 < 
0.9 or O.I2 for z > 0.9, extinction correction uncertainty Sm ext = 
0.02, low-mid 2 magnitude offset uncertainty &m l ~ m <°f f = 0.02, 
mid-high 2 magnitude offset uncertainty 5m m ~ h >°f ' = 0.04 

TABLE VIII: Forecast la constraints on the dark energy 
equation of state. (The "vanilla" parameter space is spanned 
by the six parameters (h, u m , cub, t, A, n s ).) 



the details of the CMB simulation and analysis. Cutting 
off the CMB power spectrum at a maximum multipole 



Probes 


Analysis 


Awo(2cr) 


Awi(2<t) 


SN la 


prior fl m — 0.31 


0.52 ± 0.036 


1.63 ± 0.21 


CMB + 


standard 


0.52 ± 0.027 


1.65 ± 0.25 


SN la 


fsky_~*f sk V 

in 1131191 


0.52 


1.31 




CMB / < 400 


0.50 


1.42 




1 yr. WMAP 


0.47 ±0.09 


1.76 ±0.42 




SN Sm = 0.06 


0.65 


1.75 




SN Sm = 0.062 


0.31 


0.997 




SN Sm = 


0.19 ±0.018 


0.68 ± 0.071 


CMB + 


standard 


0.20 ± 0.021 


0.37 ± 0.034 


SN la + 


WL on linear 


0.39 


0.96 


WL 


scales only 








CMB I < 400 


0.23 


0.47 




SN Sm = 0.06 


0.24 


0.48 




SN Sm = 


0.12 ±0.007 


0.29 ±0.03 



TABLE IX: Consistency checks of our dark energy con- 
straints. All constraints shown are at the 2a level, with the 
standard deviations in constraints from multiple Monte Carlo 
simulations shown when available. The "standard" analyses 
are those described in Sections IIII Bl IIII CI and IIII Dl while 
the other analyses differ from the standard ones as specified 
in the Analysis column. 



of 400, or replacing the WMAP 8 year data set with a 
simulated 1 year data set, lead to insignificant changes 
in the dark energy constraints. On the other hand, the 
constraints are very sensitive to changes in the quality 
of the supernova data. Changing the SN la systematic 
uncertainty, or its red shift dependence, leads to signifi- 
cant changes in the dark energy constraints. Meanwhile, 
the final 2a uncertainties provided by the combination of 
CMB and SN la data are nearly identical to those found 
by fixing £l m in the SN analysis alone. Thus, in some 
sense, adding the CMB data set is equivalent to fixing 
VL m as a function of wq and wi in the supernova analysis. 

We began Sec. IIV Cl bv identifying three milestones for 
comparing constraints on wo and w\ . Recall that the first 
milestone is to distinguish between a cosmological con- 
stant and a dark energy with w(l) = 0. From our fore- 
casts we find Aw(l)(2er) = 1.11, which is more than twice 
as high as the uncertainty of 0.5 needed to distinguish 
between the two dark energy models. The second mile- 
stone, to distinguish between ACDM and tracker models 
of dark energy, is also not yet reached by the combination 
of CMB and SN la data sets. We find Aw(z*)(2a) = 0.19 
at = 0.27, twice the uncertainty of 0.1 needed to reach 
this milestone. We can see from Fig. EDthat, if we shifted 
the contours to be centered on (wq = —0.9, w% = 0), then 
the cosmological constant model as well as some of the 
tracker models (wo > —0.8, w\ > 0) would be within the 
2(7 contours. Thus we cannot confidently claim that the 
combination of CMB and SN la data can rule out cither 
the ACDM model or the trackers. Finally, this combina- 
tion is far from reaching the third milestone, which would 
require that the uncertainties in wq and w\ be reduced 
by over an order of magnitude. 
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w 

FIG. 18: Forecast CMB, SN la, and WL constraints on dark 
energy. The la, 2a, and 3a contours are shown. The location 
of the fiducial model is marked by a " +" sign. 



The addition of weak lensing to the analysis resulted 
in the contours shown in Fig. 1181 with grid spacings 
A Wa = 0.06 and A Wl = 0.12. Weak lensing tightened the 
2cr constraints to Awo(2a) — 0.20 and Awi(2a) = 0.37. 
The standard deviations in the 2a uncertainties, found 
from five separate simulations, were 0.021 and 0.034 in 
Awo(2a) and Awi(2a), respectively. As above, the stan- 
dard deviations in each of the wq and wi uncertainties 
were about 10%. We checked our results using a Fisher 
matrix calculation and found the 2a constraints 0.16 and 
0.38 in wq and w\, respetively, consistent with our mini- 
mization results. Meanwhile, the minimum 2a uncertaity 
in w(z) was found to be Aw{z*){2a) = 0.097 at z* = 0.52. 

Our constraints are about the same as those forecast 
for SNAP SNe la with Planck priors, but worse by a 
factor of two than SNAP SNe la with weak lensing, as 
listed in Table IVlTTl Meanwhile, galaxy cluster measure- 
ments based on the more ambitious LSST survey claim 
an improvement by a factor of three over our forecast wq 
constraints 0. 

The analysis was repeated with two different fiducial 
models in order to assess the dependence of dark en- 
ergy constraints on the fiducial model. Moving from 
the ACDM fiducial to the fiducial model (wq = —1.3, 
w\ = 0.5) increases uncertainties in wq and w\ by about 
70% each, to 0.34 and 0.63, respectively. Moving in the 
other direction along the degeneracy curve, to the fidu- 
cial model (wo = —0.7, wi = —0.3), leads to a modest 
decrease in wq and w\ uncertainties to 0.15 and 0.28, 
respectively. Qualitatively, this is the same behavior as 
was seen with the CMB and SN la data combination. 

Further tests of the robustness of our constraints re- 
vealed that, once again, CMB data at I > 400 do not 
contribute much to the dark energy constraints (see Table 
IIX(I . Meanwhile, if the weak lensing analysis is restricted 
to linear scales, then constraints on wq and w\ weaken 
considerably, as shown in Table IIXI This is consistent 
with the findings of |39j . 

Comparing 2a constraints with and without weak lens- 




-1.1-1.05 -1 -0.95-0.9-0.85-0.8-0.75-0.7-0.65 

Wo 



FIG. 19: x 2 contours for CMB, SN la, and WL "data" sim- 
ulated using the fiducial models (»o = — 1.3, twi = 0.5) (top) 
and (u;o = —0.7, Wi = —0.3) (bottom). The la, 2a, and 
3cr contours are shown. The location of the fiducial model is 
marked by a "+" sign. 
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FIG. 20: Aw(z)(2a) as a function of z, for the data combina- 
tions CMB + SN la (upper, solid line) and CMB + SN la + 
WL (lower, dashed line). 



ing (including nonlinear scales), we see that the uncer- 
tainties in Wq shrink by a factor of 2.5, and the uncer- 
tainties on w\ shrink by a factor of 4.5. In order to un- 
derstand the contribution of weak lensing to the overall 
analysis, we compared w(z) uncertainties as functions of 
redshift for the combinations CMB + SN la and CMB + 
SN la + WL, as shown in Fig. [3U1 Without weak lensing, 
w(z) is well constrained only around z ss 0.3. Evidently, 
weak lensing adds information on w at some higher red 
shift, complementing the constraints from CMB and SN 
la. The result is a w(z) uncertainty that is not only lower, 
but much more uniform across the red shift range. 
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This improvement allows the combination of CMB, SN 
la, and WL data to reach two of the three milestones 
identified at the beginning of this section. At z = 1 
we find Aw(l)(2a) — 0.22, which will easily allow us 
to distinguish between a cosmological constant and dark 
energy models with w(l) « 0. Thus, weak lensing will 
either confirm or conclusively rule out the dark energies 
with wq + w± < favored by our analysis of current 
data. Secondly, the constraint Aw(z,)(2(t) = 0.097 at 
z* = 0.52 is tight enough that this combination of data 
can distinguish between the cosmological constant and 
tracker models of dark energy. Thus, weak lensing will 
allow us to rule out a portion of the most interesting re- 
gion of parameter space. On the other hand, our forecast 
dark energy constraints do not reach the third milestone, 
which calls for these 2a uncertainties to be reduced to 
~ 0.025. In particular, the uncertainty in w± is higher 
than this by over an order of magnitude. 

We have not included in our analysis either weak lens- 
ing tomography or galaxy cluster measuremen ts, which 
may lower statistical uncertainties even further 
These have been analyzed by others and look promising. 
They suggest that we can approach the third milestone, 
but there is no method suggested so far for pushing sub- 
stantially beyond their forecast constraints. 

Moreover, we remind the reader that the constraints 
discussed above, as well as the improvement due to 
weak lensing, are based on several optimistic assump- 
tions listed at the beginning of Sec. IIV CI Systematic 
effects, such as dust-related dimming of the supernovae, 
can increase uncertainties in wq and w\. Also, if the best- 
fit dark energy model is near the boundary of the (wq, 
wi) parameter space, then w(z) parameterization depen- 
dence will further weaken our constraints on the dark en- 
ergy parameters. In the worst-case scenario Wq + w\ « 0, 
the best fit model will be near the boundary of parameter 
space; parameterization-related uncertainties in Wq and 
w'(0) could completely swamp any improvements. 



V. CONCLUSIONS 

The purpose of this investigation has been to determine 
how well w(z) can be resolved using currently planned 
astronomical observations. The well-known challenge is 
that individual measurements do not constrain w(z) di- 
rectly, but rather some functional that depends on w[z), 
integrals of w{z), and a large set of additional parame- 
ters. 

We have shown that the uncertainties in measuring 
w and dw/dz at z = can be reduced dramatically 
by the beginning of the next decade, using a combi- 
nation of the highest-quality CMB, SN, and WL data. 
However, the remaining uncertainties, Au>o(2er) = 0.20 
and Au>i(2cr) = 0.37, will not be enough to determine 
definitively whether dark energy is inert (a cosmolog- 
ical constant) or dynamical (quintessence), unless the 
true value of w differs from —1 by significantly more 



than 0.1. Unfortunately, many quintessence models have 
|io + 1| < 0.1 and \wi\ < 0.1. 

Our numerical studies illustrate how the measurements 
combine to produce this constraint. Even the best super- 
nova measurements are degenerate under certain combi- 
nations of variations in w(0), w'(0), and fl m CMB 
constraints are degenerate along a surface in the space 
spanned by w(0), w'(0), h, and O m . However, we have 
found (Fig. [H§ that combining CMB and SN measure- 
ments effectively collapses the degeneracy in the tt m di- 
rection, leaving only the degeneracy in the u>(0), w'(0) 
plane. We have shown that the CMB contribution to 
this degeneracy breaking comes entirely from the I < 400 
region of the power spectrum. Since even the first year of 
WMAP data has reduced uncertainties in this range to 
near the cosmic variance level, we do not expect a signifi- 
cant improvement in dark energy constraints from further 
WMAP or Planck CMB data. Of course, this relies on 
our choice and number of cosmological parameters, and 
in particular, our assumption that fix = 0. If curvature 
were considered, the angular diameter distance degener- 
acy would degrade the information on Q m provided by 
the first CMB acoustic peak, and better information on 
subsequent peaks could prove valuable. Also, we have 
not considered nonlinear effects su ch a s the gravitational 
lensing of the CMB power spectra |125| , which could pro- 
vide valuable information on structure formation. 

Meanwhile, the final constraints are sensitive to the su- 
pernova systematic uncertainties. A 50% increase in 5m, 
from 0.04 to 0.06, changes the 2a constraints on w(0) 
by about 25%. After the CMB and SN data sets are 
combined, the remaining degeneracy runs along curves 
of w(z) which intersect one another near z — 0.3. The 
uncertainty in u>(0.3) is roughly 0.2 at the 2a level, de- 
pending on what functional forms for w(z) are consid- 
ered. 

To break the degeneracy between w(0) and w'(0), more 
data must be co-added that can constrain w(z) at a 
greater red shift. We have studied the weak lensing power 
spectrum as a means of breaking the w(0)-w'(0) degen- 
eracy. Our previous conclusion about the highest CMB 
multipoles remains valid; CMB data at I > 400 do not 
contribute to dark energy constraints even when weak 
lensing is considered. Furthermore, any improvements 
due to weak lensing depend crucially on our ability to 
use shear measurements on nonlinear scales. This will 
require a better understanding of systematics, such as 
intrinsic alignments of galaxies, as well as more accurate 
computations of the matter power spectrum on nonlin- 
ear scales in dynamical dark energy cosmologies. When 
CMB, SN, and WL are combined, the supernovae con- 
strain w(z) at low red shifts z « 0.3, while weak lensing 
constrains w(z) at higher red shifts (see Fig. 120(1 . leading 
to improvements of a factor of 2.5 in wo and a factor of 
4.5 in w\. 

Section HV CI began by identifying three milestones by 
which progress in dark energy constraints could be mea- 
sured. 
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• Distinguish between w(l) = —1 and w(l) « 0. 
CMB and SN alone are unable to reach this mile- 
stone, while the combination of CMB, SN, and WL 
data can distinguish between these equations of 
state at the 4.5cr level. 

• Distinguish between w = — 1 and tracker models 
with w > —0.8, dw/dz > 0. Once again, CMB and 
SN alone are unable to reach this milestone. With 
WL added, A and trackers can be distinguished at 
the 2a level. 

• Reduce 2a uncertainties in w and dw/dz to w 
0.025. Even with CMB, SN, and WL data, this 
milestone remains unreached. 

Thus the combination of CMB, SN la, and weak lensing 
is a promising tool for improving dark energy constraints. 
However, in the worst-case scenario that experiments find 
w » — 1, these three probes cannot decisively rule out 
quintessence models in which w differs from —1 by a few 
percent. 

Other works have considered highly ambitious surveys 
of the cluster abundance evolution |43j. assuming large 
numbers of observed clusters, or of weak lensing tomog- 
raphy |124| . based on multiple red shift bins and obser- 
vations of large fractions of the sky. Both probes mea- 
sure the structure growth rate as a function of red shift. 
These have the potential to improve statistical uncer- 



tainties significantly, by factors of 3-5 compared with our 
results, but c oncerns remain about systematic uncertain- 
ties jisl Il24j . Even taking the current estimated errors, 
uncertainties in wo and w\ are several percent, which still 
allows a range of plausible quintessence models. Thus, 
unless we are lucky enough to find a dark energy that is 
very different from the cosmological constant, new kinds 
of measurements or an experiment more sophisticated 
than those yet conceived will be needed in order to settle 
the dark energy issue. 
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